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1. lutroductiou 

witliiu the body are foiiatantly bombariiwl with a large repertoire of molecules 
that must be dualt with &$> potential stimuli. Most of the time, tbe;^e molecular inputs 
are measurcHl by receptors at the surface of the cells. State of these receptors are thus 
mformative un the outside world, a:id experimeutal and theuretieal biophysicists [5T] 
hav*e extensively used htformatiou theory' to estimate how' much htformatiou (m the 
Shamxou seuse can be encoded (see e.g. in a developmental context work of Gregor 
aud euwoikers [20|h 

Over longer time-scales, mfurmation proccssmg ewntually leads to a dccLsiou. that 
wt* de&ue as a cltange m physiological macroscopir behavior or in the steady states of 
a geue regulatory network. Decision is by essence a computation process, based on 
sometimes limited information. For instauee. in the well-studied example of bacterial 
cbemotaxLs, a cell might dec ide to switch behavior bet^s'oeu tumbling or swimming (S|- 
a strategy that can be wcU-explalned by a Max/Miu game theory' model It is not 
cveu always desirable to optimi 2 e mibnnation coUection if the enviroumeut changes 
too rapidly, as Illustrated by thc'iiJo taxis" strategy |23]. Other examples include 
t'ellular commitment to a given fate in response to dynamical signaling pathways [TT] . 
or decision to take action in the case of innnunc responses, characterized by binary 
Erk phospborylatiou, cytokine release, ami cell proliferation The later decisions 
are irreversible, iudicating that computation is accompamed by infurmatiou erasure 
aud thus energy dissipation (3S]. AdiUtionally, Information procH'Ssmg may be multi¬ 
tiered in order to retrieve different features of the input stlmulir rapid decision could 
discriminate betw'een ligands of different nature, while slov^'er decision could report the 
(juantity of ligands. 

In this review, \s'e will focus ou oue speciffc cellular decision where siuiilar 
conskleratious apply: the ability of T cells to discriminate very' specifically between 
self and not-self ligands. The functional significance of such ligand discruninatiou is 
(juite obvious. If a T ceU Tecoguixes'’ a (potentially single] llgaiid as foreign, a large 
set of responses is triggered to eracUcate the pathogenic Infection that generated this 
stimuli. Couversoly, if a T cell interacts with (many present) self ligamLs os self it should 
remaut quiescent to aN'old auto-lnmiune catastrophe. 

This disoiimlnation ttusk is paitlcularly daunting as T cells are constantly exposed 
to a large nimiber of molecular stimuli at once. This issue of bignalling pleiotropy 
is potout lolly a very generic problem in biology' and we will coin the term “absolute 
discriininatiou" to describe It. Midtitudes of receptors aie indeed shared in examples 
as {liffereut as BMP siguaUing, olfaction, endocrine sigualling. etc... More theoretical 
works havG suggested that orgaitizatiou of Immune repertoire requires strong overlapping 
signals (SI13S]. lu the context of tliis review', absolute dlscriininatlon thus is the specific 
aud sensltb'e recognition of foreign Ugands, independent of ligand (juantity. Absolute 
discrimlnatiou of ntassive anioimts of self vs low amount of foreign ligands is expected 
to be challenging since some self-llgands might be very close biochemically to foreign 
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Uganda. 

TLe tii'Kt part of this review will be devoted to a formal introduction to the problem 
of immune rccugultiou by T cells, priaeutiug current experluieutal umlerst ending, paat 
and present attempts to model this decision problem, and intruducmg the paradigm of 
adaptive sorting. In tbe sectmd part \s*e will introduce our eiurent model for absolute 
immune discriinination. at the cellular scale. Finally, we will discuss how tools borrowed 
from statistical physics are ueederl to understand the higher level of processing in the 
immune system, at the cellular population scale. 

2. Theoretical approaulies for absolute ligaud discriminatiou. 

T cells probe their environment in search of potential foreign peptidAi. This is <loue 
via the interaction of their T cell receptors (TC’R) with ligands (pMHC), presented 
by Antigen Presenting Cells. At a given time, these ceUs “present'’ a repertoire of 
oligopeptides {embedded within an MHC] that is representative of the current proteome 
(l.e. a mix of peptides from the self genome as well as a potential genome of the 
pathogen). The core function of T cells is to scan such repertoire anil delect the prcswjce 
of patbogeii'derived ligands ami respond, while not rcspondiiig to self'derived ligamls. 

This decision must be. by essence, al^solute in the seu.sc tliat it must be determined 
by bgaud quality (here self or nol-self}, iudepeuclently of ligand quantity (i.e. how many 
hgands there are), lu that context, deebiou has been shown to be logically all'Or>uoue, 
via binary/bistablc response in Erk phosphorylation [HSll or in NFAT translocation 
to the nucleus. A first dificulty is that there is no qualitath’e biochemical difference 
(e.g. a clenr-cut structmal tUstinction) between self-derived pMHC’ an<l pathogen- 
dcriwd pMHC, Another natural hypotltesis would be that foreign ligands lead to 
specific allosteric nio<lihcatious (conformational changes) at the level of T cells receptors, 
which \st>uld be aji ideal way to confer extreme sensitivity and specificity to immune 
recH)gnitiou. Molecular immunology has made a lot of progiess in listing all components 
imphcate«l in this early respumse, but could not find evhlence for such a direct qualitative 
sensing in the general case. Hence T cells must make a discrimination decision l>asinl on 
continuous quantitative biophysical differences between self aud not-self. explaining why 
mathematical au<l physical modeling must be called upon to tuidress how continuous 
variation in ligand characteristics gets proccfwed with alisolute discriminatiou. 

S.i, Insight from the Hfclime dogma of antigen discriT^inahon. 

Antigen discrimination in set by the hfetime of tfie aniigen/receptor complex. Tl»e exact 
molecular events associated with self/not-self ligand discrimination hy T cells remain 
elusi\t 2 . However. Ininiuuulogists, struc tural biologists aud Inophysklsts have made great 
progrt-ss to extract key parameters that physicLsLs can build upon to tackle the issue of 
specific iimuune seusing (Figure A), 

Tlie first Insight came in the 9()s svheu researidiers measured the bioph>‘sical 
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duiiadeiiatios of ligand-iec’eptor interaxliaius uKing piuiticd proteiiiK assayed uilro 
for bindliig/debiudiug (e.p. detection by Kiuface-pltusmon retionaiiee ox by oabriiuetry). 
Kersh ef cstablislicd a liler£M:rhy of bgauds with similar binding activities, where 
the life-time of the ligand-receptor complex deteiixiines tbeii ability to trigger respomse 
(Table 2 iix Pl^). Qualitatively, there exists a thresbuld of binding time (around 3-Ss) bo 
that for bgonds with a lower binding time. T cells do not respond, while for bgauds with 
bigUer binding time T cells do rwponcl (see in particular Table 2 au<l Figure 3 in ^T]). 
thus realizing absolute dlscrunination based on bindliig. Such experimentally-derived 
rule (sc^ealled “lifetime dogma" [T^) wtus well established by the turn of the miUemnmu 
such that it became the springboard for many modeling efforts. 

It sboul<l be innnedlaCely pointed out that like any dogma, this one is not absolute, 
and there are exceptions to the rule {see for a dlscmeiou of some exceptions. [TU| 
for an experimental approach on ccU pupulatioms}. For instance, it has been seen that 
some ligands could honiebow* “compensate" a small hfetime with a very high km |^> 
which ha.s been interpreted a.s an effect due to constant rebinding (101 . Mexte recent 
measmeuteuts have been carried In the context of more complete biological systems 
(e.g. T cells reading their Ugamis on the surface of antigen prA»euting cells), with 
single hgand resolution: thene brought about a correction on parameters for association 
auil dissociation rates, but cuncurred qualitatively with the previously-acqnirerl in vitro 
lueasmemeuts. 

Rcceut work by C'heng Zhu and cow*orkers La.s challenged the lifetime dogma, using a 
Molecular Force Spectroscopy to interrogate iiulividual pMHC TCR interactions, Such 
technique relics on micropipette manipulation of pMHC-coat«l beads and exquisite 
uieasmeu»eut of the force induced by the engagement with one individual TCR on 
the surface of T cells to resolve the dynamics of hgand receptor engagement. Zhu 
tt a/.'s measurements Iea4l to paradoxical results at first: strong ligand.^ that trigger 
T cells were fuuu<l to be w'eaker binder to the receptor, thus luvertiug the life-time 
dogma Subsequent studies tested how force loading on the ligaud-rcceptoi 

cuiuplex would alter its hfetime, and the mure intiutive liieraic’hy of ligauds was 
recovered, with better binders Imlutlug better signaling respoiiscs Zhu and 

colleagues thus proposed a dynamics structural model, whereby ogouLst ligands iinluce 
a conformational change in the complex (su-called catch bonds) that triggers T cell 
activation, Alternatively, non-activating ligaudh {c,g. self-dcrivetl peptide MHC) would 
not induce such confurniational change, would be rcleosetl raplilly (s^called shp bund) 
and would fail to activate a significant signaling lespouse. Hence, there would in the 
end be qualitative dlifeieuccs between activating and non-adivating ligands. This resiQt. 
while potentially establishing absolute hgand discrimination at the structural level must 
be reconcilerl with the observation that non-activating ligoucLs may be sufficient to 
activate T c'clls when proper external cues are provided in the form of cytokines (cf 
section Ultimately, the mechanical dliferences between catch bonds and slip bonds 
for pMHCVTC’R pairs, as uncovered b\' Zhu el ai. remain to bo interpreted in their 
capacity to trigger signal traiLsduction, 
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From the Physits point of view, one intriguing aspect of tliKie forco measuremuiits 
\st>ul<l he to add utetljanjctil a^peots to ligand discrimination. Understanding the 
fortres associatod with ligaud-rcceptur mteractiou.H and the coupling witli the mechanics 
of membrane doformation uxjuld be critical to account fox the diiferential potency 
of ligands to activate T cells. Such qiumtitatix'e models Lave been introduced 
l^a.Hwl on Gmzbmg-Landau equatioits coupling tin? biochemistry' of ligand-receptcn 
interactioiis with the energetic cost of membrane deformattom Such models established 
that biochemical/meLhanlcal coupling could be sufficient to physically soil ntembrane 
proteins on the T:APC' cell interface, and generate a threshold of activation. Such 
physical models generated intriguing piecUctioiis tliat were subsequently validated 
experuxientally: of note, it predicted that a family of ligands (with intermediate 
binding capacity) would abrogate the forixmtion of sc^called iixnnunologlcal syiiapse 
(a selbassembk'd bull-eye structure at the smfacc of T cells, where TC’R aggregates 
at tlie center tlie synapse, and adbesiou molecules occupy* the periphery of s>*napse). 
lu our contc*xt of immune recognition, one mmst point out that such synapse 
formation occurs dowiLstream pa.sscd the initial signaling response associated with the 
hgand discrimmation: it may constitute a reinforcing mechanism to anchor ligand 
discrimination over longer timescales, rather than the core cell-decision we are focusing 
on in this review. Recent models have explored how membrane stiffness inilueuces 
effective binding times via suptadiffu.sive effects [T]. 

Two Oilditbual lines of W'ork must be tulded to the biophysical conunflrum of 
self/not-self ligaxtd disctiinination Iw T cells. First, in the field of immunotherapy, 
researchers haw engineered T cells wdth synthetic chimeric-antigen receptors (CAR) 
whose extracellular doiualn is composed of an antibody recoguinlng a protein on the 
surface of tumors to be taigeted {e,g. CD 19 for B cell lymphoma), and whose 
intracellular domain is derived from signaling components of T ceUs {Figme [^B): 
engagement of these receptors (with non-physiologlcal Ugandn of surface antigems with 
wry laxge lifetime) ba.s been shcxwn to be necessary and sufficient to activate T t'clls. In 
fact, examples of supra^physblogical lifetimes for antigen/ret'cptor couiplexes that load 
to T ccU ac tivation were derived experimentally by in \itro evolution of the TCR/pMHC 
complex PO]. In the context of modehug early immune detection, this is relevant as 
the biophysics of ligand-receptor interactiou are wry different (with very large binding 
offiiiities), yet comsistent w'ith the lifetime dogma: antigen/receptor pairs w'ith very' 
strongly-held complexes, and wry large lifetimes are Indeerl very stimulatory'. 

Another line of experimental evidence has recently been reinforemg the lifetime 
dogma. MarkiLS Taylor fow'oikers engineered a new class of chimeric antigen 
receptors, whose extraceUular recognition unit is composed of single-stranded DNA 
[^C). Antigens for these T cells are composed of complementary single-strands of DNA 
(e.g. an ohgomer of adenosiitAs and cytokines, to as'oid secondary structures). Hence 
immune detection in that context is higliJy tunable, quantifiable and easy to ixiodeh 
it is essentially the biophysics of DNA hybrhlhation that driv'es the engagement of 
this artificial antigen receptor. In that context. Taylor et al, demonstrated that the 
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a.Hsueiatiuu rates of these artiticia] reoeptor/ligand pair were essentially constant as it Ls 
limited by the nudeatlon of double stranding between two ctmipleinentary DNA pairs. 
Howewr. the dissociation rates are highly variable and essentially dominatecl by the free 
energy of double-strand formation. Hence, these DNA-lwsed diimerio antigen ret'eptor 
and ligands recapitulate the biophysical characterfcstics of ligand-receptor mteractiuu in 
the natural immime cletection context. Most striuklngly. Taylor ct al. found that the 
hfetime dogma holds with a threshold of acti\'ation set around 3s for the lifetime of the 
antigen-receptor complex |?7g|. 

As of 2015, although the structural details of the early ewnts iu immune recognition 
by T cells remain elusive, the consensus around the lifetime dogma is thus holding and 
it is enabling physicists to iiuil<l biochemically-explicit or pheiiotypic models of good 
biological bigniheauee ^[^22]' constitutes a rich parachgm for both theoretical and 
experimental biophysical ctaLsiderations, au<l most of our dlscnssiun will be within this 
framework, 

2.9. Setting the problem for physicisU: what dots absolute immune dtscnminafio?: 
entath ? 

In recent years, quantitative inmimiology has partially chaiatleri^ed the "phenotypic 
space" of T cells aa a function of these parajueters, A "golden triangle" characterizing 
immune response can be <lrawn (TH|- Figure [^.4. The first vertex of this triangle is hgand 
speaficily, aa encapsulated in the lifetime dugma describctl in the previous section: there 
exists a:i alaolute discrimination thr^Lold on ligand binding time, around ^5 s. 

The second vertex of the triangle is hgand sensitivity. Minute amounts of ligand 
are able to trigger response. Actually, there are strong experimental evidence that one 
foreign ligand can trigger immune response (3^, bo that the physical limit of detection 
is reacherl biologically, a situation reminiscent of other famoiLs examples such as photon 
semslng Such high sensitivity might be functionally critical as the immune flj'stem 
can "snip" a pathogenic infection before it htus a chance to expand. 

The last \‘ertex of this triangle is decision speed. We know horn experiments that 
immune decision at the single cell level is taken witliin a couple of minutes . Note 
that this decLsiou time depends quite strongly on ligand concentration [2] yet it is 
relatively fast to accommodate the limite<l time T c^ells spend scanning the surface 
of one antigen presenting ceU. 

To reformulate this problem in a generic vw. imagine a cell with a given set of 
klentieal receptors is suddenly exposed to 1 Ugauds. with identical binding time r, 
We can plot in the (L, r) plane a 'tespoiise hue", characterising the bouniiary l>etween 
responding regions {"Agonibt* Uganda) and uou-respouding onc-s (“non-agonibt" Uganda). 
The life-time dogma states that below some critical time r^, fur all L. there is lio response, 
while above r,. for any L there is response. This defines a vertical line in the (L. r) plane. 
In biological terms, this is itsually called a "specific’’ response to a category of Ugantls. 
and thus we call this line “specificity line". This line corresponds to the first vertex of 
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Immuna r^cogniiion 
Miin natural TCR/pMHC 


^ Invnuna racognibon 
wlUi anlibody*basad 
chlmeflc antigar racapw 


Immune racognibon 
wiOi DNA-basad 
cblmarlc anilgari racaptor 


Figun* I: Tluee L’xaiuplos of Immuue recogiiltiou by T ch?11k. A. The natural system of 
autigcu disoriiuiimtiou by T lympliooytes relies ou the T cell Receptor, It is composed 
of an extracellular a/b domains that interact with peptide-MHC complex on the surface 
of mitigcn'prcscuting cells, and 6 intracellular <lomiuiis contaiuiiig 10 Immunorcceptor 
Tyrosme-baswl Activation Motifs {ITAMs) that get pbosphorylated upon engagement 
with uot'Self ligandSs a:id trigger T ceU activation, B. Recent developmeuts iu tbe field of 
immunotherapy iutroduccxl Cliimeric Antigen Receptor (CAR): its extracellular (lomain 
is composed of a monomeric antibody that is specific for aji antigen on tbe surface of 
the targeted tumor; its intracellular domain concatenates one ITAM^eontamhig domain 
{?.) ami adrlitional eostimulatory doutaius {CD28 and 4* IBB) to induce robust T cell 
activation upon CAR engagement with its ligand. Such CAR combines the specificity 
of an antibody-baswl recognition, with robust signaling iespuu.se. C. A DNA«based 
Chinterk Antigen Receptor was recently proposed by Ron Vale and c*ow<rtkers pS]: its 
intracellular domain is l>ased on the conoateuatiou of signaling domains -similarly to 
the C‘AR described in B, Its recognitiou platform Ls composed of a single-stranded 
oligonucleotide, that recognizes another smgle-stramled oligonucleotide by sequence 
complementarity. Rather than relying on natural pMHC ligands whose biophysical 
characteristics are not tunal)le, Vale et al.'s design cau be engineered to achieve varialile 
hfetime for the ligaml-receptor complex. Such ingenious experimental dcsigu will be 
critical to probe tbe sufficiency and limits of the lifetime dogma. 


the golden trlaugle. 

It sboiJd be notefl immediately that it is not obvious how such response cau be 
realized, especially it seems a priori impossible to haw such specific response for very 
small L (tt'e can not have response without signal!). Indeed, in any kind of biological 
settings, an obvioiLs physical limitation is that there can not be than 1 ligand 
presented at any given time. Immune cells can nevertheless trigger response when 
exposed to 1 to 3 agonist ligands; this defines a (horizontal) sensitivity line, The 
correspoucliug idealized lifetime dogma responst* line is displayed on Figme B, A thir d 
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<limeu.sioii \s*uuld be neceseary to account for tbe third clement of this triangle, epeed, 
but 1*5 not drawn here. 


2,3. Early attempts at modeling immune recognition: hneHc proo/rsadin^ 


Hi(»toiically, McKeitbau was the first to prupuse a mecliajiistie model to underly the 
early ewuts iu immune recogultiou, accoimtiug far sucL (qualitative ligand recoguition. 
The control of the fjuallty of immune respumse by a single kinetic parameter r is 
reminiscent of the famou.H Hopfteld-Ninio kinetlC'proofreadiug (KPR) paradigm, first 
proposed in the context of DNA rtplication and protem translation (HI In the 
immune context. McKeithau pointed out that subsequent to TC’R«pMHC interactions, 
the receptor does intcmally go through sowral rounds of phosphurvlatiou. C'ailing C„ 
the ligand-rec^eptur c*omplex that Las reached the degree of pLosphorylatiou in the 
easciuie. L the ligand and 77 the receptor, suuplihed cquatioms for a ctmtinnous model 
of this process are: 



- (tf. + r-‘)Co + «(L - XI C.)(J? - XI C,) 

(1) 

Cu* 

4 4 

- (0 + r-*)Cn + 0Cn-i. 1 < n < JV 

(2) 


— T + Ipc^-l 

(3) 


where is tlte association rate between ligand and receptor and o the phosphoiydation 
rate iu the cascade (here and in the foUowing, we make a meajt held approximation and 
concentrations are measued in imits of “molecule per cell’). Parameter r c*orrespoud.s 
to binding time of ligand to the rec'cptor. It appears in equations for all C„b, accounting 
for the hypothesis that after unbinding of the pMHC‘ ligand, rec(?ptor would be quickly 
dephuspLoiylated and the phosphorylation cascade would need to restart from “gromid 
wro"* {i.e. Co). 

At steady state, assuming R is far from saturation and that 0 « t~*- one can 
easily derive that the last complex C.v thus has concentration scaling as : this is 

the usual geometric dependency characteristic of kinetic proof reading with N steps. Its 
role is to amplify diflerence between ligands: with only 2 steps, for the same initial ligand 
concentration, one can get up to 6 orders of magnitude iu tlic differeijce of concentration 
of Cjv for seif hgands (r * 0.1 s ) vs agonist ligand (r * 10 s). So if immune 
decisiou Is taken via a downstream thresholding mecLanb'm, physiological concentration 
of self ligands can not trigger immime respomse^ even though one single agonist ligands 
theoreticaDy can (see response line on Figure j^D). 

However, there aie several quantitative shortcomings for a simple proofreading 
model. Fust, it is well known that to work efficiently. KPR need.s to be very 
slow, which renders it incompatible with the observed fast rt-sponse times of arlaptiw 
immune responses |2I. Second, the response line of KPR for any reasonable nmnber 
of proofreading steps with realistic decision Rp(?ed simply does not account for the 
oljserved specificity in terms of binding time r, as illustrated on Figure^ D. Finally, 
mixtures of ligand with different binding times will yield pmely additive response to a 
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Figure 2; Absolute discrimiu&tion for physicists. A: Goldou triangle for immune 
leoognitiou, B: IdeaLzed response line c’orrespoudmg to the iimuune golden triangle. 
Cl Tw<i pussible dLstributiuus for an Output variable directing immune decision, for 
uctw'orks exposed to random coucentratioue of ljgan<Ls with ulentical binding times. In 
the top example, typical t-alues for ligands with r « 3 .** and r * 10 * overlap, so that 
it is not possible to discriminate betsveeu tbese ligands. In tbe bottom example, those 
(ILstribiitious aie v^^ell separatetl so that it Ls possible to idioose a thresholding procedure 
on this Output to ensure alisolute discrmiinatiou. D: Typical response hues for various 
models discussed In this review. Decision thresltohl was adjusted to have r* 5 s, 
Parametorh for KPR and immune model from ^21' horn |30l for ailaptive sorting 
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kinetic proofreading ineeliaiiLsm, and thus would not account for more puzzling {but yet 
fuudaiuental) aspects of immune i espouse sucL as autagoniauj, where soiue nou>re«ujtis‘e 
Ugands can actually inliibit response of not self ligands [1|. 

Tbus, oue needs to augnicut the Lrailitbual KPR scheme with feedback regulation 
in order to be able to aecomit for oin “golden triangle * {epeoificity. sensitivity and speed] 
that cboracterizcs the Gorly events of T cell activation, 

2.4- /n Jiheo evolution and adaptive norting 

The moKt counter-intuitive au<l puzzbng property of innnunc response is its liigh 
specificity. The reason is that we woiQ<l e,xpect a priori that some shorter binding time 
could be “compensated” by blglrer ligand concentrations. More quantitatively, sunilar 
to the kinetic prooheading mo del, one would expect in genemJ that any output 0 of 
a general signaling would behave as O * /(i.r), where / is a mouotouic fimction of 
both 1. r. But then L«)w could we have a sharp proceffi so that, on the response curt'e. 
a sirrall decrease in r (from agonist to self) leads to a diange of L or several orders of 
irraguitude 'f 

To answer this question, we turned to tn hHuu) evolution (3S| (a review of thc' 
method can be found in [H]) , The idea is to simulate a Darwinian process on a space 
of possible models to select for absolute dLsedmination, Corisiderlng a population of 
biochemical notMt>rka (typically 30), our algorithm randomly mutates netwt>rka, PosslWe 
evolutionary moves are insphed by proofreading-based models of immune recoguition. 
and consist in addition/removal of proofreading steps, as well as additiou/iemovaJ 
of Internal phosphorylations ami dephosphoiylatkius, or of kiuase/pltosphatases. For 
selection, wc need to define a scoring or "fitness” function, and we chose to use mutual 
information as explained l>eluw. 

Assume a cell is exposed to oue type of ligands, binding time r, witli probability 
We want discrimination to be efficient over several orders of magnitude in L 
concentration. Ii» the absence of any other informatiou we diootio pf{L) to be uniform 
on a log scale, within physiological concentration range. Such a choice is also consistent 
with the well-known fact that cunceutrations in cells are distributed log-normally |3. 

Consider an Output variable 0[|]To each couple {L.r) corresponds achstributlon of 
output variable 0 characteristic of the Kignalling pathway, that we call pf{0\L }. Since 
oui problem for absolute ligand disciimination impUes immune detection independently 
of hgand concentration, let ua marginalize over all possible ligand concentrations and 
define a probability distribution for this Output associated to a binding time r : 

PAO] = J MO\L]pAL)dL (4) 

Tills probability distribution is then a pure function of the binding time of the 
hgand. Good ligand discrimination wiU be possible only if there is very Httle overlap 

t ibc udtiirv uT 0 hi luidcT sc'kt'tive preMuic in cite ulft<»rltimi aa Uiat tivulutiou I’au cUoost* tlie 
variable carryluii inaxiuiiun intbrimiilon kir diaotiuiutfttiou ms O 
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between (ILstxibutious Lorresponding to dlifeieut rs (see Figure for an example of 
dwtxibiitious witli two different ligand types). 

A practical way to use this for m silico evolution is to consider a situation where 
eellK can l>e exposed to two different types of Ugauds (say r * 3.v and r * 10a to fix 
ideas). Tu maximize) selectbu for efficient dlseriiiiktatiait. we assmue equal piobabiUty 
for observing these two rs. Then, based on computed distributioits such as the ones 
on Figure [^C, we compute mutual information between the Output and r, based on 
P^(0). A mutual information of 1 bit means that perfect discrimination Ls possible. 

This evolutionary procedure quickly converges to a very simple scheme described 
in ^ 3 . that we called adopltve soriing (FigureA). Simplified continuous equations for 
adaptive sorting are: 


Co = ^ c.){n - ^ c,) - (ox(Co) + T-‘)a,, 

Cl— Pa:(Go)Co —T *Ci, 

<P 

Co + C* 



(5) 

( 6 ) 
(7) 


Tite basis of this network is a one>step kinetic proofreading process, that is 
modulated by regulation of the phosphorylation rate 0 k from Cq to Cj. This term 
is a Alichaelis-Menten fmretiou that can be interpreted as a repression by the fir.nt 
complex in the cascade (Co) of kinase A* responsible for its phuspLorylation. Total 
contribution of phosphoiydotbn 0 k(Gd)Co therefore contains two Co dependency: a 
direct linear increasing contribution Co (substrate of the phosplioiydatbn), au<l on 
indiretl decreasing 0 k(C o) (regulation of phosphorylatiou by substrate), thus cncoiliug 
a sO'Called iucohercut fccdbrwartl loop ^S]. In the limit of high Co (au<l thiLS high £), 
those Co dependencies compensate yielding an output at steady state from equatious 

1130 


Cl * r0K(Cu)CL» 0r if Co » C (8) 

The later expression is independent from the amount of hgau<l presentctl. and then Ls 
a pure function of binding time as iUustrated on Figure]^ A. Any thresholding process 
on Cl call thiLS efficiently discriminate between binding times, B espouse hne of network 
is illustratwl on Figmej^D for a simple thr^LoIdmg proeeffi on Ci (“Adaptive sorting 
N*!”), in close agreement with the idealized lenpouse from Figmej^B. 

It should be stressed that there are two very important biochemical a.ssumptious 
related to kinase A* in adaptiw sorting: 1 , it shoulil <liffuse rapidly inside the cell and 
<j, it should iielong to a pool shared by all receptors. Sint'c any bomid receptor can 
deactivate A\ total A* thus aggregates global information over multiple bound receptors. 
A* in turns tune local state of receptors, so that the total activity of the kinase A’ is 
a decreasing function 0 a'(Co) of total Co. Overall, K effectively couples the cbfferent 
lec^eptors. wliich explains the non-linearity of the output as a fimetiou of the ligand 
concentration. 
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Figure 3; Adaptive aoiiing iictwoik.s. (A) Oue'SU’p adaptive soitiug topology evols'ed 
m ^ 3 . Output couceutmtion as a function of ligand for different vtdues of r Is 
represented. DatiUtnl blut» li»*^ Indicate thresLolding for dccieiou for deterministic 
systems used far sliuuJAtiana. (B) Adaptive sorting with more proofreading steps. 
Output coneentratiou for agonist (r * 10 .-*) mixed with ami 10^ sub tlnesholds ligands 
illustrates antagonism (response for pure ligands Is qualitatively identical to the network 
of panel A au<l thus is not showm) 


Revisiting the golden triangle« by design, adaptive sorting realizes speoificity (i.e. 
dlscriiuinates bet\s*eeL different rs) in the limit of high ligand concentration. Sensitivity 
(to small bgaml concentration) depends on the relative values of Co and C“: asymptotic 
bmit of equation [§ is realiawl only when Co ia Ihglier titajt C*. For lower value of 
Co, ipK i^ approximately constant au<1 one recovers a oue«step kinetic proofreading 
regime. So wliile discrimination very close to tbrt shold can not be fully absolute for 
small L. stronger agonists still trigger response at low concentration as indeed visilile 
on Figure [^D, Spcwl is related to low number of molecules: considering tlie immune 
example, If a T cell is able to perform detection for bgau<l c'onceutration as low as uue 
bgaml per cell, then one would expect a potentially deleterious sensitivity to stoeLastic 
ffuctuations. A natural answer to this problem is to timevaveragu 1 espouse, but then it is 
not clear any more if a quick decision can be made. In the adaptive sorting mc‘cltauisni. 
if we assmue that Ci Is activating a downstream slow output, it can be shown that 
indeed, intrimic fluetuatious can be averaged out witbin tens of seconds, which is then 
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(.'unipatible with the experimeutaJly obyervwl kminme dceb^iou time and the appropriate 
sensitivity and specificity of the TCR signalling pathway pD| . 


2,5. Antagonism 

In the previous section, wc have considered discriminatioii between ligands with diiferent 
binding times when only one typo of ligands is presented, Hoxvevcr, a more reolLstk 
immunological situation is that agonist ligands are presented simultaneously with many 
sub threshold (self) ligands. So cells not oiJy need to discruninaU* agoubsts from self 
hgands^ but should detect agonists presented within many self ligands. 

It turns out that the adaptive sotting scheme presented in Figure A does not 
perform wellr addition of few sub thresholils hgonds considerably decreases output 
concentration for the same agonist concentration. The fundamental reason is due to 
the coupling of A’ between multiple receptors dlsoaseed in previous section. Imagine 
self hgonds are arlded, calling Da conespuudiug complex (quaUty T*). then we have for 
the phospliorj'lation activity: 

and the total output concentration is 

C. + A * + -Do)(rCo + tM (10) 


For r * r,. we naturally recover the same result as before, but if t, « r, the total 
output concentratiou clearly is much lower than the asymptotic response an effect 
called antagonism, Intiutively. the "selT hgand titrates the kinase necH.’ssary for the 
proofreading step, okin to a “dog in the manger" effect dcsctibetl in other immune 
contexts If respomse is due to a thresholding effect, as a consequence, many 

more agonists hgands would be required to trigger response. In terms of respoiLse 
hue of Figme[^ B. while specificity is conseiwed, the s^'stem loses sensitivity to minute 
concentrations of ligands. 

It can be shown mathematically that antagonism is a neoessAiy consequence 
of absolute discrimination. (juaJifyiug as a '*phenatypic spamhel” ^]. Intuitivrfy. 
absolute discrimination uecessariy requires some htterual variable (simiJar to kinasu 
A’) to discriminate between ligands with different rs irrespective of their concentration, 
and as a comsequence antagonism will always occur when those internal vaiiables are 
activated by subtlm»boId ligands, A model perfuxining perfect alxaolute discrimination 
is presented in ^7\, as well as a simple categorisation of antagonistic effects, that can 
bo increased or mitigated as r —► 0 depending on the model considered. 

In particular there is a simple way to ininimixe the range of binding times with 
Htroug antagonism in the model of Figure A, by the addition of a short upstream 
proofreading cascade as discussed in Qjl]. Assuming uosv that kinase A~ is actlMited 
after r?? proofreading step.s, we have 


C», + D,n + C* 


0K(Gn» + — 


( 11 ) 
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Blit we C«, ai ran<l D„ <x r"'*'* so that Cm >> An if Ik In&h eiioiigb, even 
if luauy self Ugauds are pieseuteii. As a cuuBequence kinast* A* Ls barely iiifluenoed 
by complex Dm. A*? said before, there nevertheless still is some antagonism close 
to thieshold when r, ~ which actually constitutes a smoking gim for absolute 
(ILscrimlnatiou mechanisms (see next SH’tiou for experimental evidence and for 
a more careful study) . Antagonism is illastrated on Figure [^B for a model with 4 
proofreading steps (response line is displayed on Figure j^D). 

While antagonism is reduced l)y addition of proofreading steps, there are hom*ever 
other tiade>otfs appealing in the system; for iiLStance adding too many proofreatling 
steps might reduce the final concentration of the output too much, which creates 
sevtjral downstream problems in terms of response times verj* similar to what happens 
in McKeithan’s KPR luoilel (see Supplement of |2T] for discusnions of this effect}. 

3. Biologically realistic uiodeb for immune detection 

As we have seen iu previous sectiou, gcueric solutious to the prablem of absolute 
discrimination arc now available witlnu the sunplc adaptive sorting framework. But 
one naturally wonders if biochemical reactions iu actual T ceUs correspond to any of 
such Platonicion view. au<l. if not. if they can be related la it in any ways. 

Adaptive sorting elaborates on a small nuniber of kinetic proofreading steps. It 
should thiLS be first pointed out that many (if not all) molecular components of the 
original McKelthon model (3^ ore indeed present in actual cells. Such “realistic'’ 
model expIaiiLS its popiilaiity and its use as a blueprint for many current models of 
immune decision (3^. The phosphorylation cascade uxiuld correspond to the known 
phosphorylation of Internal Tv'rosine Ai tlvatiuu Motife (a,k,a. ITAM) eoutoiniug chains 
iu the TC’R complex [35]. R api<l dephospliorylation upon unbinding would fit tlie kinetic 
segiegatiou mechanism, specifying that generic phosphatase s are segregated only upon 
hgand-receptor interaction [I3| - 

SA. Ntyaiivt fttdhuLk and anlayumsm 

Adaptive sorting l>ased models are cxpH’ted to include an extra negative feerlback, i) 
buffering for hgand concentration to realise aiiaptation/alK»olute dLsc rimination, and 
associated with ligand ontagunlnm. Tliis has been indeed observed and characterized 
in seminal papers by Dlttcl et oi [T^, and Stefanova ei ai These papers 

establish the existence of a negative component in ligand detection via the Tyiosinu- 
protein phosphatase nou-receptor type 6 (PTPN6). also kn«nvu as Sre homology region 
2 domaimeontainiug phosphatase-1 (SHP-1), and its role iu ligand antagonism. In 
details. Dittel ci at quantified aittagoulsm by measuring T cells responding to agonists 
in coujunction with incrc^aslng eouceutratioiLS of antagonist ligands (which decreases 
the magnitude of immime response) in T cells endowed with two separate TCRs, Their 
meosmemeut also suggested that antagonism Is associated with SHP-1 association with 
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Figure 4: A siuipMed luodtO fox' immune (letectiou A: SketeL of some of tbe luulecuirtr 
plnyers aiul tbeu mterax;tiant». B: Simplified model baited ou kinetic prooCroailiag 
coinbiuod with negative feedback C:C;v coucentratloii as a fuixctiou of ligand for different 
values of r D: Total output concentration as o function of agonist (r * 10 ligand In 
presence of 10^ sub tliresbolds bgands> showing antagonUin. E; behaviour of tbe network 
when 5 total t'oucentration Ls doubled, showing collapse at bigb ligand concent ration. 
F: Experimental quantification of antagoixism. reproducer! from OVA is agonist, 
El is weak antagonist. G4 is strong antagonist. 
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TC'R aiid tl»e euaucd dephuapliurylatiou of ITAMs, Iiiipurtantly^ the antagonibtic effect 
k “iufectiona" within the ceUs: recruiLment by .SHP«1 s»preadB from receptor bound to 
autagonistic ligands even to unbound receptors, suggesting a global coupling of receptors 
via SKP-1. Involvement of SHP«1 in antagouLsm is definitely established in ^S]. They 
show in particular how moderate increase of SHP-I activity gives several orders of 
magnitude mcroaso of antagonism potency. Antagonistic ligands aro shown to 
recruit more rapidly SHP«1 compared to agonist ligands (w'bich only recruit SHP«l in a 
later time}. They finally show that agonist hgan<l specifically trigger a positive fcetlback 
loop, mediated by ERK«1. and that ERK«1 specifically inhibit SHP-I reernitment by 
the TCR, thus explaining the kinetic difierence between hgonds. 

3.9. Modeling negative feedback, approxiinatc adaptive norting 

The first model combining these diflerent aspects was pubhshed by one of its {G-A B) in 
collaboration with R.N Germain in 2005 |^. This work cumbined a kinetic proofreailing 
backbone with a SHP-l mediated feedback and an ERK-1 positive feedback. We 
inclmled most known components of the s>’stem, including different co-receptors, 
kinaift*s. and eventual phosphorylation cascade in a very complex mathematical 
framework including around 300 dynamical vaiiables. Tliis model succeeded in satisfying 
the previously described golden triauglo as a modeling target. It also estabhshes a clear 
hncar hierarchy for antaguuism, where stronger antagonists arc* ligands with binding 
times just below critical threshold r^. Most iiuportautly. the fuil-blowu model wA' 
tested with new expeiiiuents quantifying precisely antagoniam strength and decision 
time of the network, and valldatiirg predictioits from the morlcl, 

WTiile the Altan-Bounet/Germain paper could explain many experimental feat men, 
it was not clear at that time if the full complexity of the model was rwiuired to 
understand the s^'stem. Can we see the "'biological wood emerg(mg] from the molecular 
tree’* as nicely formulated by Gunawar^lena (5^? A first simplifiwl model was proposed 
in 2008 l)y Llpnlacki tt al but still contained around 40 variables, and was especially 

focasing on possible bistable properties of the system via ERK positive feedback loop. 
In 2013. we pubhshed in coUab oration with G, Voisiime, E.D, Siggia M. Vergassola 
a coiisiderably simplified Nursion of tins model [23]. focusmg on the part of the decision 
network upstream of ERK. Sciiematic of the model is displayed in Figure[7] B, Equations 
for this model are 


5 ^aCi[ST-S)-liS (12) 

s s 

(=0 (=0 

+ (i + -,S)C,-(i!> + T-')Co (13) 

C, =aCj_. + (6 + 7S)C,+,-(0 + i) + '>S + r-')Cj (U) 

C,, = ac,v-i - (6 + + T-')C« (15) 
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Like many otLei uiudelK. the backbone of thla Loaiae'gramed model relied on a 
kinetic proofn^ading backboue. Then, a simple negative feedback is included, in the 
form of the actlMitiou of a global phosplrataae called S (correspoudiiig biologically to 
SHP-I) by a single complex in the phosphorylation cascade. 

This model essentially recapitulates all observatioms made in 2005 satisfies the 
golden triangle, anil is scmi-analj-tic as detailed in the .Supplement of ^3]. In details, 
hucarity of equations to in the limit of unsaturattd receptors allow's to easily 
express conoentratiomi of all complex G’j as a fimetiou of S concentration. Defining 
> r. the roots of the chaiatteristic equation that corresponds to e<iuatiou [T^ 

Q*^ + (t + 5)r*-(0 + f) + S + v,)r (IG) 


wc can sltow that 




a7?r +1 


for 0 < j < N and : 




(l - ^-) 
kRt +1 \ / 


(17) 


(IS) 


One can then use equation that rclatee S to Cj to close the system, which gives 
a fourth order polynomial equation that can be most easily solved numerically. 

Response line of the model Ls illustrated on Figiue {“Inuumie model"). 
Interestingly, this network actually performs an approximate atlaptive sorting: it ilattems 
out the concentration of the output (7v over sowral orders of magnitude of ligand 
concentration 1. on different plateaus as a function of binding time r as shown on 
Figure 2] C , Thu.s adaptive sorting appeals to be the core piiiic iple of early immune 
detection, as could have been guessed in retrospect from iirst principles constraints 
evolved in ailico. It is known that adaptation can be performed either via feedback 
or feedforward interactions, so that the network from Figure [4] B can be seen as a 
fc'etllMick wrKion of the feedforward adapth’e sorting network presented in Figure 
Another difference with adaptive sorting as discu.ssed before is that the same klna.se 
and phosphatase is shared between all proofreading steps, while networks such as the 
ones displayed on Figure require one spociffc kinase for the step actually performing 
adaptation. Our model implements approximate adaptive sorting with a minim um set of 
sharctl uiLSpeciiie kinases and phosphatase's, and as such can be considered as optunmu 
in terms of parsimony of biochemical species used. Full stochastic simulations show that 
decision times are consistent with data from [2] and can be as low' as 30 s for strong 
agonists. Finally, only arlaptis’e sorting models with several proofreading steps such as 
the one of Figure B and Figure give antagonistic properties similar to experimental 
data Figure |4|F |23]. 

Intuititi\‘c analytical limits can be obtained w'hcu negative feedback does not 
saturate, wliich is biologically rele\iuit given the observed phenotypic variabihty (see 
pn] ami next section). For instance, if wc assume that phosphatase S (corresponding 
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to SKP-1) is in excess and tbAt ^ << 1 .in tbe limit of high Uga:id. we caji abuw tliAt 
Cl 5c 5 at \fL asing ec^ualioua [I^ [l^ and [l^ and ciiaracteristic L’quatiuu ^ gives r+ — 1 
mill r_ c< d>f5. Tbt» latter cxprc^lon of r_ is particularly iutcre;^ting, being a ratio of 
the forward phosphorylation rate 0 over tbe baeksvard depbosphorylation rale scaling 
like 5. From equation It appeals that tbe balance of these two rates essentially 
determines level of response, whicb fits tlie intuitive biological wlea that tbe main role 
of tbe uegotive feedback is to restrain progress in tlie proofreaebng cascade. 

Putting everything together, foi high concentration of agonist bgand.s. irrespectis'e 
of biiKliug times, a simple a.syuiptotic relationship holds 

C^ocl'-T (19) 

connecting ligand roucentratiuu L and Cjv. with N the number of proofreading steps, 
So if iV is high euuiigh, C,v decreases with bgand concentration L in this model. This 
seems rather coimter-intuitiw, and predicts that if feedback is strong enough, immune 
response at high agonist ligand concentration could thus disappear (sec Figure]^ E, full 
stochastic treatment con be found in (131) Indeed, we tested and v'erified this prediction 
in cells with high level of SHP*1 |23|. Furthermore, if SHP-1 level is increased abu\‘e a 
couple a few-fold, iiegath'e feedback essentially dominates for all ligand concent ration, 
auil response is fully abolished, a “digital effect* indeed experimentally observed in [IH] 
and first predicted wltli tlic AJtan-Bonnet Germain model 

.9.5, Beyond absolute diacrimination ? 

Other modelling strategics for immune recognition by single cells have been proposed. 
Major differeuces in must models is that r depeudeucy is not considerc*»I os sltarp os wbat 
is assunii*!! here, and mudulatiuu by other parameters is considered. For instance it has 
bevu Kiiggested that there might actually be an optimal time of interaction (i.e. respoiLse 
woulil decrease at higher rs) ^]. This effect is debated, since strong binding TC’R.s 
evolved m wire, with rs of the order of 50 to 100 s conversely yield an exceptionally 
potent respoiiBO [TS]. Lever et al recently reviewed au<l compared 5 faniiUes of models 
[T2| . including the one from previous section, and huggestetl that a mid-range r can 
be simply explainfti by a kinetic proofreading with "limited* signalling’*, meaning that 
the sigualling complex in the cascade gets inactivated with constant rate, However, 
this model does nut include the negative feetlback descrilx'tl here, and thus do not 
display antagonism. Another issue of iuttrest is the regime of fa.sl "on rates"’; for such 
hgands, response is experimentally observed even for low rs. It has beeu suggested that 
competion bet weeu biuding/rebiudiug and diffusion of receptors on the membraite would 
actually give an effective higher r (23' which therefore wmld not ncc'essarily contradict 
the life-time dogma but rather cumpl^sxify it. A recent review [T?3 summarkes muleculai 
players of the flj'stems and lelationships to proofreading mechanisms, discussing in 
particular the role of co-receptors au<l various kinases, and possible other mechanisms 
that could explain better high sensitivity such as serial triggering by one hgaud. 
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Other theoretical uiodcls Law explored biopliyHical Ibrrita or more reliiied 
(.“omputationB. For iiutance, a relatetl (luestiou is the niiniiauio tlreoretieaJ decision 
time for such process ( “decision on the fly"). As said in introduction, a trado>off 
between accurax;y and precision is expected, and tills is a potentially acute issue in 
an iuiuiiine context. If declsiun take's too much time, a foreign ligand might be gone 
before proper response is activated. If accuracy is decreased, there is potent id for autO' 
immune response. Exact results for this problem have been recently obtained b>' .Siggia 
and Vergassola ]!^. They study the problem of detection of change of composition 
of ligand mixtmes. It tmns out that Wald sequential probability test ratio on the log 
likelihood of a sequence of binding e\t?nts can be u.swl to take decision. Strikingly, simple 
phosphorylation networks reminiscent of network controlling early immune detection can 
naturally implement biological versions of this test |^, An important aspect to optimize 
decisiou time is that decision here Is maile ‘dyuamicaJly ’. while networks performing 
decision liere essentially u'ork at steady state (see nevertheless Supplement of 1301 fur 
<lis( Uivsioas of transient behaviours). Such dynamical sensing couhl also Itelp in the 
regime of very small ligand conceutrationti. maybe in conjunction with serial triggering 
[TQ| . An even more general problem Is cliemodetectiou of Ugaiids of different qualities 
in fluctuating environments, where mixture composition can vary wit time. If there 
is a high probability of observing ligamls just subtlircsbohLs. antogoiusm is helpful to 
essentially buffer Bpurioii.s flucluationa of ligands |3T]. reminlst'ent of what is observed 
in olfaction Finally, some models Lave recx’ntly explored theareticalJy the possible 
network topologies ma.'dmizing information transmission on ligand mixtures, anil haxG 
pointed out the need of both proofreading and internal feedbacks [THl 157] , 

4. Ligand recognition by a population of lymphocytes: more is different? 

4-1. Tatkliuy cell-tO’cell muabilily in immune lesponses: -slatislical physwH to tfie 
rescue? 

There aie still many challenges to fully understand early immmie respou.se. In particiJai. 
once a eeU has been activated, it appear.s that further procAjslng occurs at the immune 
population le\t?l. Modeling the early events in immune detection potentially is very' 
topical for statistical physics as it involves the accounting of ccll-to-cell variabilities, 
modeling immime ruipuuses with distribution of cells, and the testing of their fimctionaJ 
signiheauce. 

Starting with the seminal work of M. Elowitz cl at [T7| a subfleld of biological 
physics has grown to address the emergence of cell-to-cell varjability||} Extrinsic 
variability stems from many sources (e.g. all source's that do not relate to the 
stocliasticity and luw-copy number in biochemical reactions); they can include 

uoce that. iJibauf^h pbysiotKU arc toad lo dc<soribc Kuch vuriAbOltY as 'noise'', bd^td ou (lieir 
n'[)rcseiitatiuu ii&ug m LouKCNiu oquatiou. (Ills ti^rui rcuMius cou/usiUK lor uiodt biolugl&Ci^ b(X'oii«e 
of uegativu (cjiiaotatioa: wu lo use the more neutral Uriu uf c’ell-iu-oc*U vuriabiiiiy 
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heteiogein?ity in the cellular euviroiuijeut, as well as epigenetic vaiiation and inetabolic 
fluetuatbib> from eeU'tocelL Intrinsic cuntiibutions can l>e best studied in an ieugenic 
population of ccUs and are the physical cunse<{ueurcs of stochastidty in the biochemical 
reactions. As a cousct^uence. ceUs can display broad distributions of phenotypes based 
on the heterogeneity of abundance of key regulatory proteins (receptors ^ kiiiaHcs and 
phosphatases in the context of signal traitsductiou network; transcription factors in the 
context of gene regulation). 

First forays to adihess the functional significance of t‘ell-tt>cell liability in 
biological systems were focused on bacterial response. From chemotaxis {the abihty 
to orient motions in gradients of nutrients or chemokiiies) to competence the ability to 
acquire uew genomic materials), researchers demonstrated that, indeed, varierl levels 
of key proteins could map into varied phenotypes ^TJIS]. Such observations were used 
at first as uew quantitative constraints to vahdate biochemical models of biological 
regulation. 

Concomitantly, these observations demoustratotl that a “mean-field' measurement 
aud model of a population of cells might have serious short coming when predicting 
global 1 espouses. One example where such c'ell-to-cell variability was foimd to be critical 
Ls in the study of bacterial antibiotic resistance. Fur example, Balaban and coworkers 
introduced a microllukhc device to track the proliferation and death of bacteria under 
antibiotic treatment ^ [S3]. A sub population of isogenic bacteria were foim<l to resist 
to aiitibioticB, simply by being a different metabolic state compared to their sister cells 
at the time of exposure to antibiotics. Such process of distributed response based on 
distributed phenotype at stimulation time was duhcrlbetl as an optimal strategy to tune 
respoiLses in a fluetuatlng enviromnent. by matching probability of phenotypic switcliing 
to probabihty of envirumucntal changes [13]. Indi*«I. there is such a fimdamental 
mismatch betVkUeu the liccessary response time (ceUs must exclude antibiotics on very' 
short timescales) and the evolutionary constraints (it \s'ould take a large amount of time 
aud a very low probabihty for cells to generate a solution to the problem of antibiotic's 
resistanc'e). that cells are better off diversifying tbeu phenotype pre-emptively such that 
a solution is readily accessible when the antibiotic perturbation applies. 

^.5, Phenotypic vanab'iliiy of T cel/ ligand dwcr:Tn:»mb'cn, 

Similar constraints are at play in the context of the imiuune response. There Is a 
similar disconnect betv»'eeu the dynamics of biological problems at stake (eradicating a 
fast-replicating fast evolving pathogen vs. generating an adaptive immune response). 
In particiilai. recent mcasiuements by the Jeiikins X: Davis lab |5TJE§] have beautifully 
illustrated the uumlier of constraints for a good immune response. Lymphocytes can 
rapidly proliferate (by factors of 10^ to 10^) and relax back to Iom' numbers for the 
memory' pool. Such explosive proliferation is critical to match the challenge posed by 
fast-proliferating pathogens. However, surprisingly, the nnuiber of T cell clones that 
can recognize a specific antigen (e.g. ilu peptide) Is very small, with iO to if)D clones 
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pur individua] (mouse or liuman}. This Is partlmlarly relevant as pathogeus caji be 
“cunning'’ an<l attempting to invoke multiple esc'ape inechaniKins. Hence ^ rather than 
a deterministic adaptation of the immune response to the pathogenic challenge and 
optimization of detection at the cellular scale, one ca:i conjecture that the immune 
system relies on some degree of statistical randomness in order to diversify rc>spouses 
within on Isogenic population of lymphocytes. 

Tlte functional mnlerpinnings of such ceU-tocell variability can be readily detected 
in the context of tbe early ewnts of Immime detection. Indeed, researchers have been 
using the va»t panoply of antibodies to quantify pxotein and phuspL(>protein levels 
in cells, as well as the single-cell resolution of cj'tometry (using fluorescence-baswl or 
luast^spoctroiuetry-bahcil approaches), in order to quantify the ceU-to-cell vuriabihty of 
responses to external stimuli. In a nutshell. If the abundance of protein X Is liniitlng 
in the activation of Y into Y* (e.g. phosphovviatiou). then measuring and correlating 
X Ncith Y* abunrlauces at the suiglc cell level will reveal heterogeneity of the response, 
Such cell-t^cell variability analysis has been carriwl out lu the context of immune 
detection to demoiLStrate the sensitivity of ligand disci iiuinatlou to varied levels of 
signaling componenth (e.g, C'D8 and .SHP-1) (5^, as well a.s the sensitivity of T cells 
response to the cytokine lL-2 to varied lev'els of cytokine receptors (e.g, C’D25, C‘D122, 
auil CD132 a.k.a. IL'2Rn. lL-2R;y and 7r) [11|< Note tliat such parameter sensitivities 
were first predicted from the dynamical model of thu signaling ca.scades at play, and CVll- 
to-Cell Variability Aitalysis (CCVA, a new methodology introduced in [I^) validated 
these prwiictious quant itatlv*ely. Sunllarly. single-cell analysis has been carried out 
in the space of phospho-pruteuis by the Pc'er & Nolan labs to quantify the strength 
of ct)unections within tbe TCIR signaling pathway ^]. There, analysis of single-cell 
measmeiiteuts using overall, lesolutlon of immime iespuu.se}4 at tbe inrbvidual ceU level 
ha.*i highlighted the large phenotypic variability in tlic signaling response of individual 
lymphocytes: such observations must then be interpreted functioually to map out how 
T cells diversify their response to optiiuize its detection capabilities. 


^.5, Dtnviny I'tiialle immune rtspoiUieH from uncehaftfe Jcaponse.H of indimdual tefis. 


The cell-tChceU variability of lymphocyte response* presents a challenge fur our current 
c‘ell-centric understanding of self/not-self dh>crinnnation in the innuune system. Indeed, 
if each individual T cell makes a rapid. sLaip yet utterly variable decision to respond to 
a hgand (see , one could autkipate many "mistakes" whereliy T cells would respond 
to self tissues and trigger an auto immune disorder. 

Wc conjecture that it is the integration of the respou.ses< of individual T cells over 
longer tiniescalea (> hour) that may conect for thu ^sloppinuss" of Individual T cell 
on short timescales (as minutes). Such integration can be carrier! out tlirough cell-cell 
conmmnicatiou.H e.g. through secreliou ami consumption of evtokiues. In particular, 
K. Tkmh ei ol. measiued expcrimeutslly bow much cytokine accumulates in the 
Hupematant of T cells that were activated m Tfiiju pS^* ^ surprising scaling law for 
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this /£ —2„uu oapiieity with the regards to the quantity of antigena 4^{pMHC) and the 
siie of the T cell popiiJatiou is observetl 

(#(pA/HC))”-“.(.Vr„„.)-'’‘ (20) 


Hence, T cells are able to output cytokiitcs In a uear«Unear scalable manner across 
four decades of #(pA//fC): this is quite a striking ol)ser\atlon when considering the 
limited dynamic range of Imlividua] T cell respam* and more generally of any biological 
system. Moieowr, the independence of output with the number of T cells Nraiia in the 
system is remarkable con.siderii\g that cytokine secretion is a priortan exteirsive variable: 
this finding implies tlrat the owrall capacity of the system (in terms of maximal cy'tokine 
concentration) Ls an inteirsive variable. Deriving an intensive output (t,e. an output that 
is independent from the size of tlie sj'stem) fur a population of T cells has been observed 
in very relalerl experiments Hart ci measured the cell expansion of a population 
of 004"^ T cells after actbution through then antigen receptor pathway Hart ei ai/s 
model f«>cuB«l on the cytokine IL»2 whose divergent fimctlou (cell proliferatiou and 
apoptosis) would explain such robust system-size-independent output. 

Understanding mcclranlstlcalJy the emergence of such scaling laws Implied revisiting 
our classical biodjemi<:a) imderstauclirig of the IL'2 cytokine pathway. It was well 
established that T cells would essentiaUy shut down their IL-2 productiou as s»>on 
an they secreted and built up a pool of shared cytoklue. Biochemically, this Implied 
that the IL-2 output of a population of T cells would have a cellirig of 10 pM. as the 
CHrncentration of cytokine required to induce IL-2 sigiiAliiig and sulrsequent shutdown 
to IL-2 secretion ^3* However, detailed quantitation of the biocliembstry of the system 
unraveled a negati\t 2 cross-talk from the antigen signaling response to the IL'2 response 
pathway ^3- measuring the phosphorylation of the STATh transcription factor 

dcrwrrstieam of IL-2 sensing. Tkwh et at imcovcrcd a surprising convolution between 
resporrse to pMHC antigens au<l /L-2 cytokine: 


#(p5r.4r5) sc 


i^{IL^2Riihc] 


1 + 


K, 






( 21 ) 


where #(-Y) represents the number of X within the e'en. ^(/L-2i?ci) is the uiunber of a 
chain of the IL — 2 receptor, au<l #(/L-2I?^/7c) is the number of complexed 
and 7c receptors on the surface of T cells hr a nutshell. ^(/L-2i?r») gets expressed 
in activated T cells to complete the sigualing #(/i-2i?i^/7r} pair into a signaling 
complex, upon 71-2 binding [T2]. Overall, this equation encapsulates the rouv*olution 
betwwn local and global responses (i.c, pAfHC abundances ajjd 71-2 concentration 
respectively): this convolution Ls particulaily significanl a.s it regidates the off awitch 
for 7L-2 production. This in important because pSTATh in turns regulate IL-2. 
which meaiLS that there is a feedback between the local stimulation (by pA77fC ) au<l the 
global readout (by 7L-2) through the regulation of STATh phosphorylation, Owrall, a 
conibination of an incoherent feedforward loop was sU«iwn to be necAjsary au< 1 sulfic ient 
to explain the scahng law in the accirmulatiou of 71-2 in the miheu. 
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Siidi initia] attempt by Tkauh tt ul whs l>a.se(l uu dotaiJoil bbcbeuiicaJ 

modeling with explicit biadiemistry being implemented. Additional experiments 
rev'ealwl additional functional relevTiuco for the secretion of IL-2. as a global legulatut 
of self/uot'Self discrimination; VoLsiime ei at probchl the proliferation response «)f a 
population of T ccEs containing multiple clones of diverse specificity {?.e. a polyclonal 
population of cells), TTiLh study demonstrated that strongly activated T cell clones could 
induce tbe activation an<l proliferation of neighboring weakly activated T cell clones . 

this previously midacmncnted pLenomenon of propagating T cell activation was termed 
'T cell co-optation". Similarly to previous studk's (H an experimentally-derived 

blocbeniically-explicit model of such multi scale iittegratiou was introduced to account 
for such l>inphocyte cooptation. Overall, it was demoustrated that the adihtion of 
autigeu signals (read locally thi'ough the TCIR pathway) ajid cytokine signals (read 
globally through the IL-2 pathway) decides T ceU fate Such observations expand 
ov'cr longer timescales (days) an<l larger timeacaJe» (lymph nodes) what iitdividual T 
cell can contribute in terms of immune response. This is particularly significant as it 
demoustrated that T cell activation cau no lougor be cousldered as a local property 
of autigeu recognition, but gets decided by integrating multiple cues. The practical 
implication should not be understated when manipulating the response of a polyclonal 
population of T ceUs to force v^'eak T cell clones to respond to antigens (e.j. In the 
context of tinnora) may hold the key to tumor eradication. Future efforts wiU require 
phenomenological coarse-graining to allow better imclcrstanding of immimc recognition 
at the level of the s>’stem. 

Recent technical dowlopments to monitor T cell responses at the indisbdual cell 
level are enabling researchers to track the early events of immune activation, one cell at 
a time. Monitoring tbe (bfferentiation of indivicinal lymphocytes wiD certainly acc'elwate 
our study of the immune Hj'stem, in particular when stochastic effects and phenotypic 
Nuriability are necessary to explain the divt'rsification in immune detection. Yet, in the 
context of the study of the innumw system as a whole, it is the c'oUective properties 
of the ceEs based ou their cytokine counnunicatioiLS and competition for antigen that 
shape the overall immune tespou.se. 

Finally, ou the pure theory side, some recent questions inspired by statistical phj-sics 
include Low the innnmie Ugand landscape Ls matched by leceptors/c'ellular diver.Hity. 
auil how collective decisioua can be made within this framework, without entering into 
specific signalling details, A simplified piobabilistic model of T-coU activation has lM?en 
derived using extreme values IheuriAs (modelling binding euergs' of TC*R.s to single 
hgand). from wliich prolmbUity of false positive can be computed[7|. and coUectis’e 
decision svas shown to improve witli the help of ceU to cell cummunicatious, Mayer ci at 
\t 2 ry recently proposed a general framework, predicting that more receptors ate needed 
for rate antigens, and strong cross-reactivity to CHirrectly cover the full antigen landscape 
pB| . The future wiU oertainly involve more statistical-physics based approanhes that 
tackle the large number of lymphocytes and focus ou such emergent properties. 
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Figure 5: lumiune rocuguition over long tlmescaJca may involve collective 
syudirouizatiou of T cell activation, \Vlijle T ccUs perform an absolute (immune) 
dLscrimlnatiuu at the individual ceU level, tbeii heterogenous exprc-ssiou of sigualing 
(.H)nipoLeuts leads to unreliability in their thresLohl of activation. Moreover, 
experimental svork on lymphocyte slguallug demomitraterl that iudividual T ceUs 
respond in an all'Or-none (digital) such that their output contains limited information, 
Additional meclianiKnm must be at play to “proofread” individual immune recognition 
and/or generate plastic inunune responses. Feedback legulatioms based on cytokine 
secretion and consumption constitute an efficient solution to the Umitatious of individual 
T cells (nUESHSOilSTI. In particular. T cells can rely on such cell-to-cell commuuicatious 
to achieve higher levels of immmie recognition, B. Tkat’h et ai found a surprising scaling 
law. whereby the maximiun concentration {[fL — 2]„uu) of the lL-2 cytokine released 
by a population of T cells scales almost linearly with the ajiiouut of antigeus that is 
present in the st’stem, practically independently of the nuinbt-r of T cells present in the 
system. C. Su<’h analog bcaling at the population Ibvel was foim<l to derive through a 
coherent feed-forward loop of Type 4, usiug the uomenclatiuc introduced by U. Alon 
oud cowurkers [45] , 


5. Wrapping up. 

To conclude, we presented an ch'erview of some phj'sics-inspircd studies in immmiology, 
focused on early recognition by T-Cells. Specifically, we discussed how the core problem 
of the adaptive immune response (the recognition of pathogeu-detlved antigeus) must 
be adrlressed with quantitative models. Indeed, experimental results have now well 
establislicd the biophysical underpinnings of self/nol-self discrimination namely the fact 
that small increases in the lifetime of ligand-receptor complexes learl to large increases 
in ligand potency, Additioually, cells can respond very sensitively and with great speed, 
Quantitatively reconciling these three experimental aspects (schcalled golden triangle) 
is a theoretical ehalleuge with great conceptual and practical relevance (e.g. when fine- 
timing T c'L’U activation is rcHjuhed, a.s In cancer immunotherapies). We reviewed the 
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curtfiijt states of theoretical models, l>uil<bng upon the original kinetic proofreading 
scheme. In silico evolution of biochemical networks satisfying the golden triangle 
ujjtavelecl a muiuual model that can reconcile all aspects of ligand disc^riniination. 
In particular, a proximal negati\'e feedback (associated with the activation of a 
pJjosphatase) was found to be critically relevant to abrogate responses to self ligands 
(ewn in large quantities) while allowing responses to not-self ligands (even in small 
quantities). Tlie functional pay-off of these models is to "predict" the existence of 
antagonism iit immune recognition as well. Finally, we discussed liow these models 
of ligand dLscrimination by T cells create new challenges in terms of understanding 
the phenotypic variability of Isogenic populations of T celb, m in terms of acr^iunting 
for the quantitative response to antigens when measuring T cell acti\'ation over long 
timescales. Similar collaborations between experimentalists an<l theoretical physicists 
will remain fruitful to expand our quantitative understanding of T cell activation to 
more complex issues in immunology (role of regulatory T cells, tuning of responsiveness 
according to infiainmatory milieu etc.). More generally, we hope that these fmtdamental 
issues of immunology will spark the interest of statistical physicists, as the derivation 
and manipulation of large-scale immune response from the local activation of incUvldual 
T cells remains poorly understood at the theoretical level. 
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